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We perform a detailed NLO analysis of the combined CCFR xF% and F2 structure 
functions data and extract the value of a s , parameters of distributions and higher- twist 



(HT) terms using a direct solution of the DGLAP equation. The value of a s (Mz) = 
0.1222 ± 0.0048(exp)±0.0026(theor) is obtained. Our result has a larger central value 
<3\ ■ and errors than the original one of the CCFR collaboration due to model independent 
parametrization of the HT contributions. The dependence of HT contributions on the 
QCD renormalization scale is studied. 

Oh! 



In the recent years interest to the problem of the extraction of the high-twist-terms from 
the analysis of different deep-inelastic scattering (DIS) data was renewed, mainly due to 
the possibility to model these terms in different processes using the infrared-renormalon 
(IRR) technique (see e.g. Refs. [0-0 and, especially, Ref. Q for the review). 

On the other hand, the experimentalists improve their data precision and achieve, 
sometimes, a percent level of accuracy. For example, very precise data on xF% and F2 
from the vN DIS experiment, performed at Tevatron by the CCFR collaboration, recently 
appeared 0J§]. The CCFR data on xF 3 were analyzed in Ref. |J in the next-to-leading- 
order (NLO), and with an approximate next-to-next-to- leading order (NNLO) corrections. 
For the latter the NNLO QCD corrections to the coefficient function |nj were taken into 
account. The NNLO corrections to the anomalous dimensions of a limited set of even 



non-singlet moments |TT| were also taken into account. The NNLO corrections to the 
anomalous dimensions of odd moments, which are not still explicitly calculated, were 
obtained using smooth interpolation procedure proposed in Ref. [12| and improved in 



Ref. |]]. The aim of Ref. || was to attempt the first NNLO determination of a s (Mz) 
from DIS and to extract the HT terms from the data on xF% within the framework of 
the IRR-model |3[]. Alongside, the model-independent extraction of the HT terms was 
made, similarly to the analysis of the combined SLAC-BCDMS data |13| , which was 



performed in the NLO approximation. Theoretical uncertainties of the analysis of Ref. 
II were further estimated in Refs. [14] in the N 3 LO approximation using the method of 
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Pade approximants. It has been found in Refs. PJTil that the inclusion of the NNLO 



corrections leads to the decrease of the HT terms, so that at the NNLO its x-shape 
variation is closer to zero. 

In this work we completed the fits of Ref.[|15[], performed the NLO analysis of the CCFR 



data on the structure functions F 2 and F 3 with the help of a QCD DGLAP evolution code, 
developed in Ref . [TO . It should be stressed that the code |TB[ was tested using the proce- 



dure proposed in Ref. |T7j and demonstrated the accuracy at the level of 0(0.1%) in the 



kinematic region covered by the analyzed data. Our fits were made in the NLO approxi- 
mation within the modified-minimal-subtraction (MS) factorization and renormalization 
schemes. The effective number of flavours rif was chosen to be rif = 4 for Q 2 less than the 
definite scale Mf and increased to rif = 5 at larger values of Q 2 keeping the continuity 
of a s [[0|]. The value of the effective matching scale M 5 was varied from M 5 = m;, to 
M 5 = 6.5m&. The latter choice was advocated in Ref. on the basis of the DIS sum 
rules consideration. The dependence of the fit results on the choice of the matching point 
gives one of the sources of theoretical uncertainties inherent to our analysis. 

The leading twist terms xF 3 LT (x, Q) and F 2 LT (x, Q) were obtained by direct integration 
of the DGLAP equations for non-singlet, pure-singlet, and gluon distributions that were 
subsequently convoluted with the coefficient functions. In order to provide the straight- 
forward way for comparison of our results with Ref. [|16[], the initial reference scale for 
pQCD evolution Ql = 9 GeV 2 was taken. The boundary conditions at this reference scale 
were chosen in the form analogous to the ones, used in Refs. ||||: 



xq 



NS 



(x,Q )=r} NS x bNS (l-x) CNS (l + 1 x)- 



Ans 

for non-singlet distribution, 
xq PS (x,Q )= Vs x bs (l-x) Cs /A s 
for pure-singlet distribution, and 

bo i 



xG(x,Q )=VGX bG (l-x) CG /A G 

for gluon distribution, where A^s, As, and Aq were defined from the partons' num- 
ber/momentum conservation and other parameters were fitted. 

The expression for the xF 3 and F 2 that includes the HT contribution looks as follows: 

xF^x, Q) = xF 3 LT ' TMC (x, Q) + F 2 HT (x, Q) = F 2 LT > TMC (x, Q) + 

where F 2 L 3 T,TMC (x, Q) are F 23 (x, Q) with the target mass correction applied. We used 
the model independent HT-expression, i.e. H 2>3 (x) were parametrized at x = 0.,0.2, 
0.4, 0.6, 0.8 with linear interpolation between these points. 

The account of the point-to-point correlations of the data due to systematic errors 
can be crucial for the estimation of total experimental errors of the parameters (see, in 



particular, Ref. ||18fl , where the value a s (Mz) = 0.1180±0.0017 (stat+syst) was obtained 
as a result of the combined fit to the SLAC-BCDMS data with HT included). The 
systematic errors were taken into account analogously to the works of Refs. Hl6| , [T8|| . The 
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Figure 1. The high-twist contribution to the structure functions F 2 F 3 . Full circles 
correspond to the fit with renormalization scale parameter fcfj = 1, empty circles - to the 
fit with ku = 1/4, squares - to the fit with k^ = 4. 



total number of the independent systematic errors sources for the analyzed data is 18 
and all of them were convoluted into a general correlation matrix, which was used for 
the construction of the minimized \ 2 . In addition to the point-to-point correlation of the 
data due to systematic errors, the statistical correlations between F 2 and xF% were also 
taken into account. The account of systematic errors leads to a significant increase of 
the experimental uncertainties of the HT parameters and the shift of their central values. 
However, even in this case, there is a definite agreement with the results on HT-behaviour 
of Ref. ||, obtained in NLO. Moreover, these results do not contradict the IRR-model 
prediction of Ref. ||, since fitting IRR model parameter A 2 to the data, we obtained 
A' 2 = -0.10 ±0.09 GeV 2 , compatible with A 2 « -0.20 GeV 2 and A 2 « -0.1 GeV 2 , given 
in Ref. § and Ref.§. 

Performing the trial fits we got convinced that the introduction of the factor (1 + 72) 
into the reference expressions for the the gluon and singlet distributions does not improve 
the quality of the fit and does not change the value of a s . Also, we fixed parameters 
Ins, bs and be at zero because this increase the value of x 2 by f ew units only while 
X 2 /NDP remained less than unity. 

The results of the fit on H 2 (x) and H^{x) parameters are given in Fig. 1. One can notice 
that, comparing with the fit to xF 3 data only from Ref. ||, the HT parameters errors are 
decreasing. Within the errors, the parameters that describe the boundary distributions 
are compatible with ones of Ref. ||. The Hs(x) coefficients are in agreement with the 
NLO results of Ref. |J and the behaviour of H 2 (x) qualitatively reproduce the HT 
contribution to F 2 that was obtained from the combined fits to the SLAC-BCDMS data 
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When the matching scale M5 was changed from nib to 6.5m&, the value of a s [Mz) 
shifted down by 0.0052 and, hence, the theoretical error in a s (Mz) due to uncertainty of 
b-quark threshold can be estimated as 0.0026. This uncertainty is in agreement with the 
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results of the fits to the CCFR data obtained within the so-called spline MS prescription 
pl|j with the help of the Jacobi polynomial metod |22| . One more source of the theoretical 



uncertainty due to the truncation of higher QCD orders was evaluated following the way, 
which was proposed in Ref. In accordance with their procedure, one can introduce 

renormalization scale k R into QCD evolution equations in the way, given below for non- 
singlet evolution: 



dxq a s (kjiQ) 
dliaQ n J x { 

. ® s (k R Q) 



2tt 



P q N q S ' (1 \z) + (3 P q N q s ^(z) \n(k R )} }- z q NS (x/z, Q) 



where p NS <(ty and P NS <W denote the LO and the NLO parts of the non-singlet splitting 
function. The dependence of the results on k R would signal an incomplete account of 
the perturbation series. The shift of a s (Mz) resulting from the variation of k R from 1/4 
to 4 turned out to be only 0.0007. At the same time one can observe a simultaneous 
variation of H2${x) (see Fig. 1). This effect can denote, that, in fact, the fitted values 
of if 2,3 (3?) incorporate the higher order QCD contributions (NNLO and beyond, c.f. Ref. 
[1)0)- ft i s interesting, that the same effect was also observed in the analysis of charged 
leptons data [^3|, i.e. it cannot be attributed to a specific feature of the CCFR data. 
The interplay between these contributions and genuine power corrections does not allow 
for their unambiguous separation. The value of ot s is strongly correlated with the fitted 
high twist, that leads to increase of the a s error and, consequently, in our fit with simul- 
taneous determination of a s and high twist contributions the theoretical error in a s due 
to truncation of higher QCD orders is merged into the total experimental errors. Having 
taken Qq = 20 GeV 2 as an initial scale, we checked that our NLO results were quite stable 
to the variation of the initial scale. However, we do not know what will happen at the 
NNLO, where the stability to was observed only starting from Ql fa 20 GeV 2 [14]. 



The final value of a s in NLO with the account of theoretical uncertainties is given as 
a s (M z ) = 0.1222±0.0048 (stat + syst)±0.0026 (theor.) It differs a bit from the NLO value 
a s (M z ) = 0.119 ±0.002 (stat+syst)±0.004 (theory) obtained in the CCFR analysis [§. 
The increase of the experimental error is due to that CCFR group used model-dependent 
form of the HT contributions, while we considered them as the additional free parameters 
and extracted them from the fit. 

It should be stressed, that the scale-dependence uncertainty of the NLO results dras- 
tically minimized to the value of 0.0007 after taking into account HT corrections. The 
decrease of this uncertainty was also recently found in the analysis without HT terms 
after taking into account 3-loop splitting function [2~4~|. Thus we think that more careful 



NNLO DGLAP analysis of the DIS data with HT effects included is now on the agenda. 
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